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The generalization of Density Matrix Renormalization Group (DMRG) approach as implemented 
in quantum chemistry, to the case of non-orthogonal orbitals is carefully analyzed. This gener- 
alization is attractive from the physical point of view since it allows a better localization of the 
orbitals. The possible implementation difficulties and drawbacks are estimated. General formulae 
for hamiltonian matrix elements useful in DMRG calculations are given. 



I. INTRODUCTION 



The DMRG method, first introduced by White 000 m solid state physics, has been very successful in calculations 
of low dimensional quantum lattice systems, and has become recently an interesting tool in quantum chemistry 
0, IS E EJ EH ■ The method is attractive because of the good scaling properties of the calculation time with the 
size of the system . 

In the present work we first list the main steps of DMRG algorithm as implemented in quantum chemistry. All previous 
implementations were restricted to the orthonormalized molecular orbitals (MOs). While rendering all algorithms 
much simpler, the standard MOs may be not quite appropriate from the physical point of view, especially if we try to 
follow the original DMRG ideas that assume the separation of the system into well defined physical blocks. From this 
point of view the use of bond structures, e.g. valence bond orbitals [1 1| can be much more attractive as it allows for 
clear interpretation of resulting wave functions. Even if it is not clear that this can improve the results numerically, 
we consider below the theoretical aspects of the use of non-orthogonal orbitals as applied to the DMRG algorithms. 
The article is organized as follows: first we describe essential steps of the classical, orthogonal, DMRG algorithm. In 
section ITTT1 we describe the generalization to the case of non-orthogonal orbitals. In appendices we give the full list of 
relevant formulae. 



II. CLASSICAL, ORTHOGONAL, DMRG ALGORITHM 

We do not describe here the details of the DMRG approach, it can be found in many reviews 0j 0, • We also do 
not describe the details of implementation of DMRG in quantum chemistry, see 0, Q ■ Rather we concern only with 
practical, computational, aspect of the algorithm. 
The main points of DMRG procedure are: 

• The parametrization of the wave function 

• List of 'primitive' matrix elements stored 

• Adding orbitals 

• Calculation of Hamiltonian matrix elements 

• Diagonalization of Hamiltonian matrix 
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• Expanded space states construction 



A. General form of the wave function in DMRG 



In DMRG we divide the whole system into A and U subsystems (two subsets of the full set of MOs in quantum 
chemistry). The wave function is represented as 

v = j2 c v\ A i®uj> (i) 
u 

where \A] > describes some state of subsystem A and similar for U. The tensor product Aj ® Uj is automatically 
antysymmetrized due to the presence of the creation operators (see 8] for details) . The conservation rules imply for 
nonzero C'u 

n a {A I ) + n a {U J ) = n a (2) 

where n a is fixed by the total number of electrons and total spin projection, the same being true for (3. If there 
present additional spatial symmetries, the corresponding conservation rules must be also satisfied. 
One important drawback of the above expression is that the wave function given by Eq.Q is not necessarily the 
eigenfunction of total spin operator S 2 . We actually never considered the application of corresponding projection 
operator to get the spin-pure DMRG wavefunction, maybe it is not really as difficult as the S 2 is the combination of 
creation/annihilation operators whose matrix elements are in general available in DMRG (see below). 



B. State pair types and matrix elements to be stored 



To find the wave function given by Eq.([Q) we must diagonalize Hamiltonian matrix in the space span by \Aj ® Uj > 
state vectors. To build this matrix all we need are some matrix elements for pure A and U blocks that we call 
'primitive' matrix elements. The set of these matrix elements depends on the relations between number of a and (3 
electrons. In the following when we consider particular matrix element 

< IJ\H\KL > (3) 

indices / and K refer to A subsystem while J and L - to U subsystem. 

Let's denote with the symbol a creation operator referring to block A and with a creation operator referring to 
block U. The electronic part of quantum chemistry Hamiltonian can be written as 

H = H A + Hu + H AU , (4) 
where Ha, Hu contain only a, or u, operators and the interaction term is given by 

Hau = \^ x y\Pl) ~ ( a; 9|ro)][ui T u y Ta| )T a gT + u^u^a^a,;] + (^^[u^u^a^a^ + u^u^a^a^], (5) 

pqxy 

where p, q indices refer to orbitals of A block, x, y to U. 

Then, given I and K, the Table [I] lists the 7 cases (others can be obtained using < IJ\H\KL >=< KL\H\IJ >) for 
which the matrix elements Eq.© are nonzero. Please note that the JL-case is inverse relative to that of IK, i.e. the 
case of LJ is the same as IK. 

With this definition, the 'primitive' matrix elements needed to build the full Hamiltonian matrix, depend on the 
'case' given above. The complete list of these matrix elements is given in Appendix^ together with resulting matrix 
elements of the total Hamiltonian (see below). Here we just note that all ME's have the typical form of 

< I\a pa ...\K> (6) 

with one to three creation/annihilation operators, or the linear combination of such ME's. 
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TABLE I: 7 cases for /^-relations 



Case 



P 



nf = n% + 1 
nf = Uk - 1 

n? = n£ - 2 

n? = < 
1/ =n%-l 



B S 
n I = K 

< = <- 
= rC - 



nf = nj - 1 



C. Calculation of Hamiltonian matrix elements 



After we have defined the stored 'primitive' matrix elements in the previous subsection, the matrix elements of 
Hamiltonian matrix look like 

< IJ\H\KL >= J2 fi{IK) 9i (JL) (7) 

i 

where type and range of index i depend on the IK -type. The full set of matrix elements is given in Appendix ^ 
where the bi-orthogonal orbitals are assumed, see next sections. For the standard orthogonal case the dual orbitals 
coincide with original orbitals so the tildes in the formulas of Appendix^ can be simply omitted. 



D. Adding orbitals 



Adding of orbital(s) is the most important and basic step of the whole DMRG procedure. Given the set of A states 
(the same is valid for U, so we do not speak of it further), we add few orbitals to A and build the new bigger set 
of states including all possible occupations of added orbital pattern. Indices Ii, Jj, etc. refer to the new enlarged 
subsystems. The subsequent step is the recalculation of the 'primitive' matrix elements in this new set of states. The 
only practical difficulty here is that the number of cases to treat grows considerably. 

In our early implementations we recalculated the 'primitive' matrix elements after adding orbitals and stored them on 
disk. This led to significant reductions in performance. Therefore later, to save disk space and improve performance, 
we changed our algorithms as not to calculate and to store the extended set of 'primitive' matrix elements but rather 
to build the Hamiltonian matrix in the extended space explicitly using the non-extended 'primitive' matrix elements 
and some additional molecular integrals. Also, the highly optimized -BLA5'|l2|, \l4 . Il5l llq routines are used as 
much as possible to make the code faster. Note that this part of code can be easily parallelized. The full list of 
relevant formulas is given in Appendix [BJ The tildes can be safely omitted again there. 



E. Hamiltonian diagonalization in the orthogonal case 



Once the Hamiltonian matrix elements in | Ai®Uj > basis are available, the next task is to diagonalize this Hamiltonian 
matrix. We can not apply the explicit techniques like Jacobi diagonalization, because of quite big, tens or hundreds 
of thousands states, dimension of typical problem. Therefore we arc forced to use the direct technique, which is based 
on the ability to calculate the so called a- vector, i.e. action of Hamiltonian 

y = Hx (8) 

given vector x without explicit storage of Hamiltonian matrix H a p. In the usual formulation with orthogonal orbitals, 
Davidson method is the most efficient. As well known, it consists of an iterative procedure that starts from some 
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reasonable guess and then expands the set of trial vectors up to say 30-50 that provide good approximation to exact 
eigenvector. The procedure can be seen as a minimization of the Raylcigh quotient 

E=^^, (9) 
< x|x > 

in which one applies the Newton method and replace D = diag H in the Hessian matrix. 



F. New states construction 



After diagonalizing Hamiltonian matrix in the extended set of vectors, the next and the last task for the whole DMRG 
procedure is the construction of reduced subset of states describing A and U in the extended orbital set. This is done 
by finding the best possible approximation to the calculated eigcnstate of Hamiltonian matrix when the number of 
states retained is fixed to a number (M) smaller than the total dimension spanned by \Ii ® Jj >; this approximation 
is obtained by minimizing 



In practice it is done by performing the Singular Value Decomposition (SVD) of the matrix Cjj (see Eq.(2.1)) which 
gives the components of the wave function in the Ai ®Uj basis, and then selecting the largest singular values and 
the corresponding states. 

After we know the expansion of selected M states in terms of \Ii > for A and | Jj > for U, we must recalculate and 
store the 'primitive' matrix elements for these new states. Again it is done 'on the fly', expanding the matrix elements 
to newly added orbitals and contracting at the same time with coefficients produced by SVD procedure. For reference 
purposes in Appendix [C] we list all expanded matrix elements. These formulas are closely related to that of Appendix 

m 



III. NON-ORTHOGONAL ORBITALS 



In this section we consider how the algorithms of the previous section can be extended to the case of non-orthogonal 
molecular orbitals. 

We start by describing general techniques to work with non-orthogonal orbitals. In principle there exist a general 
Lowdin formula that expresses the matrix elements of one- and two-electron operators between two determinants 
build with non-orthogonal orbitals as the sum of appropriate matrix determinants. It is this necessity to evaluate 
many matrix determinants that renders this approach completely unusable in practice. The practical way to work 
with non-orthogonal orbitals consists in introducing the so called bi-orthogonal (dual) orbitals 0, 0] . 



A. General remarks on the eigenvalue problem with non-orthogonal bases 



First of all, we cannot use the variational approach for the simple reason that we are unable to evaluate directly the 
diagonal energy expression which we demonstrate in the following section |20j. 

In order to compute the variational expression (Eq.(2.9)) for the energy, we must be able to express the vector |\& > 
in the biorthogonal orbitals; if we were able to find this expansion, then the computation follows more or less the 
standard root with insignificant changes since the one-electron and two-electron integrals loose their transformation 
properties. However in practice to find and use such expansion is an almost impossible task. In fact when the orbitals 
of A and U do overlap, a given state of the block A would become a linear combination (of nearly full CI size) of 
different A and U states. Therefore the whole strategy of DMRG becomes unapplicable. Still we will show that some 



5 



modern non-variational approach can be used with more or less success. 

Let us introduce some notations that will make clearer the following discussion. Given a set {fa} of non orthogonal 
orbitals, the dual or biorthogonal basis fa verifies the condition: 

<fa\fa>=S ij (11) 

Let S denote the overlap matrix < fa\<f>j > ; given two one-particle states \a >= J2i al \4>i >■> \b >— J2ib l \fa >j we 
can define two different scalar products: the "physical" scalar product 

<a|6>=^aVSy (12) 

ij 

which will be denoted by Dirac brackets, and a simple "numerical" scalar product defined by: 

(a\b) = (13) 

i 

which will be denoted by round brackets. Furthermore for any state \a >= ^ i a l \fa > we define the dual state 
\a >= J2i o^ifa >! therefore we have the identity: 

<a\b>={a\b) (14) 

The concept of dual state can be generalized: a dual many particle state is the state built with the same CI (or 
DMRG) expansion coefficients but with {fa} orbitals replaced by dual ones pi) . 

A nice property of the biorthogonal formalism is that the "dual" annihilation operators defined by 

aia = ^^aicyi^S^ 1 )^ (15) 

k 

obey the usual anticommutation relations |17| : 

{a£5jv}+ = £ ffT fiy (16) 

The one-particle and two-particle integrals are computed with the fa wave functions on the right and the dual fa on 
the left. For instance: 

(ij\kl) = / fa(l)fa(2) — fa(l)fa(2)d ri dr 2 (17) 
J r 12 

In non-variational approach instead of looking for the minimum of energy in the trial subspace we project the energy 
eigenvalue equation onto the 'dual' subspace. To find the current approximation to energy and vector given the trial 
vector set we solve the following non-symmetric eigenvalue problem 

HijCj = EC, (18) 

where HTy =< Xj|.ff|xj >= (xi\H\xj); because of biorthogonality < Xj|xj >= (xj|xj) = Sij. If expansion set becomes 
reasonably full for the description of target state, this non-variation approach must produce reasonably good result, 
which is also demonstrated by some experience with non-orthogonal Full CI calculations. 

In principle, if we want to set up a variational calculations with non orthogonal orbitals, we need the knowledge of 
the transformation operator T from normal to dual vectors : 

x = f x (19) 

Then the variational energy becomes 

E = <x|g|x> = < x|f T ff|x > = (x|f T H|x) 
<x|x> <x|f T |x> (x|f T |x) 

An iterative minimization procedure for this energy expression can be easily written down; however, it is almost 
impossible to follow the variational approach. Even the simple square < x|x > of the physical norm of a many 
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particle state |x > which appears in the denominator of Eq.(3.10), is a huge expansion of determinants, since the 
scalar product of two simple configurations gives rise to a determinant; neglecting the spin indexes, and denoting by 
|0 > the vacuum, we have: 

< 0\a in ....a il af 1 ...afJO >= det(< (j> ih \(j> jk >) = det(S lhJk ) (21) 
As we shall see in the next paragraph, the same kind of difficulties exist for DMRG calculations. 



Once the non symmetric Hamiltonian matrix Hij is constructed, we are left with the problem of using Davidson 
method or a similar direct meyhod. In paragraph 2.5 we have seen that the Davidson method can be considered as 
a minimization procedure of the expectation value (2.9). This procedure looses any meaning for a non symmetric 
matrix, since the minimum of (2.9) does not coincides with the lowest eigenvalue. However the Davidson method can 
be seen from a different angle, and can be used for solving the eigenvalue problem of non symmetric matrices. Of 
course the orthonormalization can be performed only with the "numerical" scalar product (3.4). 



As well known for non symmetric matrices, there are other good methods that can be tried. The simplest one 
is perhaps the two-sided or non-symmetric Lanczos method that treats on an equal footing the Krylov sequences 
{x, Hx, H 2 x...} and {x, H T x, (iJ T ) 2 x...} generated by the matrix H and its transposed H T . Let us denote by V and 
W the two supspaces generated by the two Krylov sequences. It is easy to find two biorthogonal bases {vj} and {wi} 
in V and W, in such a way that the matrix (wi\Hvj) is tridiagonal [l9j . 



More complicated methods, like the two-sided Jacobi-Davidson method, which makes use of biorthogonal basis like 
the two-sided Lanczos, can also be considered 0] . The use of orthogonal bases is generally more numerically stable, 
while two-sided methods (and biorthogonal basis) may provide faster convergence. 



B. Variational calculations with non-orthogonal orbitals in DMRG? 



Here we show that variational calculations are completely impractical for non-orthogonal orbitals. We consider the 
overlap matrix just. The Hamiltonian matrix would be even much more complicated. So the problem is to evaluate 

< IJ\KL > (22) 



Actually the problem is not just to evaluate Eo Q22JI but do it efficiently. For that we must factorize this expression to 
something containing only pair quantities, as it is anyway impossible to store and use M quantities, corresponding 
to Ea (|22fl . This means that we would like to have something like 

< IJ\KL >= x l {IJ)x 2 {KL) + x 3 (IK)x 4 {JL) + x 5 (IL)x e (JK) (23) 



It is rather easy to show that such factorization is impossible in general case. We show that by considering the most 
simple example when all states /, J, K, L are determinants, and the number of electrons in / and K states is the 
same (we consider say a electrons, (3 being similar). Note that while for orthogonal orbitals the block structure is 
automatic, i.e. number of electrons in / and K (and also in J and L) MUST be the same to get non-zero overlap, it is 
NOT the case for non-orthogonal orbitals. Then the overlap, according to Lowdin formula is given by the determinant 

< IJ\KL >= det ( SlK SlL j (24) 

\ SjK Sjl J 

where S'-matrices are buildee by the corresponding orbital overlaps. Assuming the general case when matrix Sjk is 
non-degenerate and using the Gauss formula we have 

< IJ\KL >= det(S IK )det(S JL - S JK S^S IL ) (25) 



which has a complicated mixed structure and cannot be represented in the factorized form Ea l|23|l . 
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C. Non- variational calculations 



As we have seen that the only way is to perform non- variational calculations, the rest is more or less straightforward. 
To evaluate the matrix elements of Hamiltonian we use the same expressions as for orthogonal orbitals with minor 
modifications due to loss of integral symmetry. The full list of matrix elements to build the Hamiltonian matrix if 
given in Appendices lAl and iBl Having Hamiltonian matrix elements (or better corresponding er-routines), we solve 
non-symmetric eigenvalue problem as described before. The only remaining point is the reduction of states. We 
cannot calculate the reduced density for the same reason for which we cannot calculate theoverlap matrices. Actually, 
the concept of reduced density matrix is not clear in the case of non-orthogonal A and U. The matrix 

M 

Pu> = CijCpj (26) 
j=i 

is not the reduced density matrix of one block and does not verifies Tr[p] = 1 . 

However we can follow a simpler approach: we can abandon the idea of approximating the wave function using the 
physical norm < x\x > and we can use the following norm induced by the simple expression (3.4): 

\x\dmrg =< x l x >=< x l^ x >= (x|x) = N 2 ( 27 ) 

i 

Of course the norm defined by Ea. (|27|l is a true norm, so it is positive, it is zero only when x is zero ecc. Then we 
can find the 'best' approximation using this norm, which as before leads to SVD decomposition of the wave function 
coefficient matrix C/j. Again when the orbital overlap matrix resembles unit, this should give results close to those 
obtained using the true physical norm. However the minimization with the "numerical" norm can give good results 
in more general cases. 

Let us consider in more detail the practical side of the calculations. Let us define OP(~<-> no ~) the operator that 
interchanges normal dual orbitals in all the one- and two-electron integrals. Clearly we have: |22j | 

OP(~^ no ~)[< IJ\H\KL >] =< KL\H\IJ > (28) 



This means that all 'primitive' matrix elements not containing integrals are the same for IK and KI. Therefore in 
practice we can keep the same number of cases (7) and just calculate two copies (one normal and one dual) of those 
primitive MEs that contain integrals, i.e. Hjk for case 1, and f x for cases 3 and 4. Then we can use Eg . 1)281) to 
calculate inverse cases. The attention we must pay is that when we have LJ instead of JL referring to primitive 
MEs, we must take the dual copies. Considering that in practice the most of storage is taken by f pq primitive matrix 
elements of the case 1, which remains unchanged, the bi-orthogonal case has essentially the same memory- and-disk 
requirements as the standard orthogonal procedure. 



APPENDIX A: NON-EXTENDED DMRG MATRIX ELEMENTS 



We start by explaining some notations we use in the below tables. The phases as as follows: the A is created first, 
and U next. Within the state the a- (spin- up) electrons are created first, and j3- (spin-down) next. So the 1 1 J > state 
where / refers to A and J to U, can be written as 

\u >= |/(T)/(|) J(T) J(l) > (Ai) 

As to letters, / and K will refer to A block and J and L - to U . Orbital labels p, q,r,s, . . . will refer to A orbitals 
while x, y, z, ... - to U. The states are ordered as follows in the program: J < K if: 
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• rij < n K 

• rij = n^ K and n" < 

• rij = n K and nf — and symm(I) < symm(K) 

Here we list the matrix elements needed to be kept and the corresponding contributions to the matrix elements of 
total Hamiltonian. As noted before we consider only the case I < K. The opposite case can be obtained from the 
modified 'Hermitian' property (3.19). For the specific case 1, when J > L in the below expressions, we can use the 
property: 

f:f(JL) = OP(~~ no ~)[/;/(£J)] (A2) 

Table |n] gives all matrix elements. 



TABLE II: Hamiltonian matrix elements for DMRG 



Case (I < K) 


MEs to keep 


< LJ\H\KL > 


J <L1 


Comments 


1 


< I\H A \K > 
f£ q (IK) =< I\h+ a a qa \K > 
f^ q (IK) - similar a <-+ f3 


Sjl < I\H A \K > + 
8 ik < J\Hu\L > + 
52 wy9x * y (IK)f2 y (JL)+ 


YN 






9 - (J/0 = E P9 fZ(IK) [(Ml*!/) - (pvI&j)] 


J2 xy g? y (IK)fg y (JL) 




Not kept 




+ E m /A( JJ£ ')(mI 5 ») 






(evaluated 




g^ y {IK) - similar a <-> /3 






when needed) 


2 


/pg(JJT) =< J|a+ a a ?/3 |/sT > 

s. y (JiiO = -E WI /«( J *')(p«lw) 


T, xy 9*v(IK)U(LJ) 


N 


Not kept 


3 


fp(IK) =< /ISpal/f > 

/^(J/f) = E p , r < I\a+ p a P i3a ra \K > {ir\qp) + 
E, ;p>r < I\a+ a a pa a ra \K > [(ir|gp) - (xp\qr)] 
" 9*{1K) = f x {IK) + J2 p fp{IK)h§ x 


(-l) n K + n K + 1 { 

J2 p f p (IK)f p (LJ)+ 
J2 x g x (IK)f x (LJ)} 


N 


Not kept 


4 


like 3 with a «-» /3 




N 




5 


f pq (IK) =< / a pa a,a -K" > 
gxy(IK) = J2 q<P fpi( IK ) \{xp\yq) - (yp\xq)} 




N 


p > q only 
x > y only; Not kept 


6 


like 5 with a <-> /3 




N 




7 


fp q (IK) =< I\a pa a ql 3\K > 
g*y(IK) = j: pq f Pq (IK)(xp\yq) 


as for Case 2 


N 


Not kept 



APPENDIX B: EXTENDED DMRG HAMILTONIAN MATRIX ELEMENTS 



In this Appendix we give the matrix elements of Hamiltonian in the extended basis, i.e. when we have added the 
orbitals to all four states /, J, K, L. The states over those added orbitals are simple determinants, which we denote as 
patterns i,j,k,l. The number of cases thus grows enormously if we count for all possible occupations combinations 
for states /, J, K, L and patterns i, j, k, I. The fact that states over i, j, k, I are complete, i.e. are simple determinants, 
allows us to reformulate the problem by shifting everything to A block, i.e. we call the 'combined' pattern i + j as 
new i and k + / as new k. The only thing that remains after we have diagonalized Hamiltonian matrix is to cast 
this combined representation back to original split one. This is just the phase and indexing problem that is done 
very efficiently. Thus we have now orbitals added only to A block, as i and k patterns and thus the number of cases 
becomes reasonable. The phase convention is as follows: 

\iu >= |/(TKUWTMI) J(T) J(i) > (Bi) 
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The orbital labels a, b, c, . . . will refer to orbitals of patterns i and k. In below expressions the sums over those indices 
a, b, c, . . . is assumed when not given explicitly. Table UTTI lists all the matrix elements. For the explanation of quantities 
/ and g see Table ITT1 The 'new' case IiKk means the relations between number of electrons in li and Kk blocks, in 
the same way as before. The overscore on the case index mean inverted case, i.e. li > Kk. 



APPENDIX C: EXTENDED 'PRIMITIVE' MATRIX ELEMENTS 



Here we give the working expressions for 'primitive' matrix elements in the extended basis li, Kk. The notations are 
the same as in Appendix [B] As cited in the text, these matrix elements are not stored. Instead they are calculated 
'on the fly' and then contracted with vectors coefficients produced by SVD procedure to give final matrix elements 
for the reduced M state space. The full list of nonzero matrix elements is given in the Table ITvl 
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TABLE III: Extended Hamiltonian matrix elements for DMRG 



Old case (IK) 


New case (IiKk) 


(Jz) < (Kk)? 


J < L? 


< HJ\H\KkL > 






VN 


YN 

I IN 


S lk < IJ\H\KL > +5ik5 jl < i\H\k > + 

5 la kJlK T,ab9ab( JL )fab( i k P) + 

&% l3 k, 3 &IK J2ab 9ab(JL) fab(iak a ) + 
S i0 k 5jL J2ab9ab( IK )fab(' i o l k a ) 




1 


2 


Y 


N 


5i K {-l) n * J2nh < ia\&t a \k a >< ifj\a bfj \k fi > g ab {LJ) 


1 


3 


Y 


N 


{Sik [5 ipkp (E a 9«(LJ) < i a \a aa \k a > + 












J2x f*( LJ ) Ea,6>c < ialataabaacalka > [(ab\xc) ~ {ac\xb)] \ + 










J2 x f*( LJ )12ab c < *c,|a C a|fea >< i/3 \a+ g a b/} \k/3 > (ab\xc)_ 


+ 










S ii3k l3 Y,cx9xc(IK)f x (LJ) <i a \a ca \k a >} 




1 


4 


Y 


N 


(-l)< +1 {5 IK [8 laka (E a 9a(LJ) < ip\a a p\kp > + 
E* f*( LJ ) T, a , b >c < i^atpabpa^kp > [(ab\xc) - (ac\xb) 
E*/*( iJ )Ea&c < h\o-c0\k/3 >< i a \at a a ba \k a > (ab\xc) 
$i a k a J2 cx 9L(IK)f x (LJ) < i \a C 0\k >} 


+ 


1 


5 


Y 


N 


SlKSi (J k [} Ea>6 < ict\a a0 i&ba\kct > 9ab(LJ) 


1 


6 


Y 


N 


8iK8 ta k a J2 a >b < i P \&a(ia b j3 fc/3 > g ab (LJ) 


1 


7 


Y 


N 


Si K (-l) n * J2„b < i^aa^kc X if]\a bfj \k fj > g ab {LJ) 


2 


1 


YN 


YN 


S.J L (-l) n * +1 Ysab < *«|Oaa|fca >< ifil^kp > g ab {IK) 


Z 


o 
z 


V 
I 


1ST 
IN 


SikJ2 xy 9xy(IK)f xy (LJ) 


2 


3 


N 


Y 


*iafca(-l)" S E c < i?\*t \h > E x U(JL)g xc (IK) 


2 


4 


Y 


N 


V? (-!)"* E c < ic\a ca \k a > J2 x MLJ)gc X (IK) 


3 


1 


YN 


YN 


{5 JL (E a 3a(/K) < i a \a+ a \k a > + 
J2 p fp( IK )J2a, b >c < ic,\at a a£ a a ac ,\kc, > [{ba\cp) - {ca\bp)]\ + 










J2pfp( IK )J2 ab c < i<*\at a \k a >< iplatp&a^kf} > (ba\cp) 


+ 










S hk E cp 9c P (JL)f P (IK) < i a \a+ a \k a >} 




Q 
O 


9 


NT 

IN 


Y 


5i a k a (-ir* +n " + < +1 < ip\a%\kp > E cp f P (IK)g P c(JL) 


3 


3 


Y 


N 


(-l) n k+ n k+ n K+ n K+ 1 { 

Sik hr p f p (IK)f p (LJ) + J2 x g x (IK)f x (LJ)\ + 












6 icka Eab < -i^lpO-b^kfJ > J2 px fp ( 1 K ) fx ( L J ) ( £ P\ 2& ) + 










5 h k fi Eat < ialaiccabdka > J2 px fp ( IK )f* ( LJ ) [( £ p\^ b ) - (ap|x6)] 1 


3 


3 


N 


Y 


E a>b < i<*\a+ a a+ a \k a > J2 px f p (IK)f x (JL)[(ap\bx) - (bp\ 


ax)] 


3 


4 


Y 


N 


(-1)<+"k+< J2 ab <i a \a+ a \k a > 
< ifj\a bf j\k,3 > J2 px f P (IK)f x (LJ)(xb\ap) 




3 


4 


N 


Y 


(-l)< +n K+< +1 j: ab <i a \a+ a \k a > 
< ip\&tp\kp > Epx f P (IK)f x (JL)(bx\ap) 




3 


5 


Y 


N 


<5 i/3fc/3 (-l) n ^ + < +1 E c < i a \a ca \k a > J2 p f P (IK)g cp (LJ) 


3 


7 


Y 


N 


s la kA~i)^ +n * + <j: c < i ^0\h>i: p fp( IK )9pc(LJ) 
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Table 3 continued 



Old case (IK) 


New case (HKk) 


(Ji) < (Kk)l 


J < L? 


< HJ\H\KkL > 


4 


1 


YN 


YN 


(_l)»S+»£+»5r+i (£ 3a(/tf) < i/S |o+ Ifc^s > + 
J2 P fp( IK )J2a,b> c < iplatpatpaa^kp > [(ba\cp) - (ca\bp)]j + 










J2pfp( IK )J2 a bc <i f ) \ d -tf3\ k f) >< i a \a£ a a aa \k a > (ba\cp) 


+ 










Si a k a J2c P 9c P (J L )fp(I K ) < hp\a c p\ip > j 




4 


2 


Y 


N 


<5 i(3fc ,(-l)"- + < +1 E c < iA\k a > J2 p f p (IK)g C p(LJ) 


4 


3 


Y 


N 


(-l)<+-K+<+^ ab<i/} \a^ \ k0 > 
< i a \a ba \k a > J2 px f P (IK)fx(LJ)(xb\ap) 




4 


3 


N 


Y 


(-!)<+"%+< <i0 \a+ g \k > 
< i a \a+Jk a > J2 PX f P (IK)f x (JL)(bx\ap) 




4 


4 


Y 


N 


(_l)»fc+"fc+™K+ n K + 1 { 

Sik (j2 p fp(IK)fp(LJ) + J2 x g x (IK)f x (LJ)\ + 












& if> k n J2at < ic*\a+ a a ba \k a > J2 px f P (IK)f x (LJ)(xp\ab)+ 










5 ick a J2ab < ipKpVbrikp > J2 px fp(IK)f x {LJ)[(xp\ab) - (ap\xb)]j 


4 


4 


N 


Y 


E a >b < ipfit&pfo > J2 px fp( IK )MJL)[(&P\bx) - (b P \ 


ax)] 


4 


6 


Y 


AT 
IN 


6 ia kA-lT« +n « +n * +1 E c <^Mh > E P fp(IK)g cp (LJ) 


4 


7 


Y 


N 


<W-l) n ^+< +1 £ c < ^\a ca \k a > ^ p f p (IK)g cp (LJ) 


5 


1 


YN 


YN 


SjL5 i/3 k J2a>b < *a OfcaOia fca > g a b(IK) 


5 


3 


Y 


N 


S igke (-l)<+< +1 Y,r. < ^\a+ a \k a > f x (LJ)g cx (IK) 


r 



r 




Y 


ivr 
IN 


Sik ZL>„ 9x V {IK)f xy (LJ) 


6 


1 


YN 


YN 


5.JLSr a k a J2a>b < i^tA^P > 9ab(IK) 


6 


4 


Y 


N 


5i a k a {-l)< +1 Y. c < i ^U k P>T JX U L J)9cMK) 


6 


6 


Y 


N 


Sik J2 x>y 9xy{IK)f xy (LJ) 


7 


1 


YN 


YN 


Sj L (-l) n * +1 £ ai) < i a \a+ a \k a >< ip\a+ B \k s > g ab (IK) 


7 


O 


Y 


N 


*ia*a(-i)" s E c < h\&tg\kf} > E x MLJ)g*c(iK) 


7 


4 


Y 


N 


5i B k g (-irZ + < +1 E, < i a \&L\k a > Y,J*{LJ)g cx {IK) 


7 


7 


Y 


N 


8ik J2 xy 9xy(IK)f xy (LJ) 



13 



TABLE IV: Extended 'primitive' matrix elements for DMRG 



Old case (IK) 


New case (IiKk) 


(Ji) < (Kk)l 


New MEs 


1 


1 


YN 


IjAa c ttA I r rja , 

s i a k a ga b (IK) < i/3|a+3O 6 0|fc0 > ] 

f a,l3 (HKk) — S i. f a,l3 (IK) 
f" b (HKk) = S IK 8 ifj k l3 < i a \ai a a ba \k a > 
fPjIiKk) = 5i K 5 ia k a < i/3 la^abslka > 


1 


2 


Y 


fab(HKk) = ( — l)™* < ialajLlfca >< | O-b/3 | > 


1 


3 


Y 


fa(HKk) — 8 IK (—l) n 'i< + ' n K 8 iaka < i a \a aa \k a > 
fJHKk) = fcjf/Jftliy o" (/K)<5 iflfcfl / H (i a A; a )l 


1 


4 


Y 


fa(HKk) = (-l) n « + < + <(5 ict fe cv < i/3|o ^|fc/3 > 
fJHKk) = (-l) n £r+ n £r [5 Jir /x(ifc) + (-1)"* T gL(IK)Si a k a fa(i0ks)] 


1 


5 


Y 


fab{HKk) = SlKSipkpfabiiaka) 


1 


6 


Y 


fab(HKk) = 5lK5i a k a fab(i(ikf3) 


1 


7 


Y 


fab(HKk) = (-l) n <=<5/if < i a \a aa \k a >< i/3\a b g\k0 > 


2 


1 


YN 


Hf^kk = (-l) n ^ + 1 T, ab 9ab{IK) < i a \a aa \k a >< i fj \a+ g \k fj > 


2 


2 


Y 


f pq (HKk)=5 ik f pq (IK) 


2 


3 


N 


f x (KkIi) = 5 iak J-l) n K +n K+< J2 c g xc (IK) < if,\&2g\kf> > 


2 


4 


Y 


f x (HKk) =6 i0k0 (-l) n ^ +nP K +1 Y Jc g cx {IK) < i a \a ca \k a > 


3 


1 


YN 


Hf° Kk = (-ir^ + < +1 [E p f P (IK)f p (ki) +E a 9a(IK)f a (ki)] 
f c JHKk) a = (-l) n K +n K + 1 8 iaka < i a \a+\k a > fJIK) 


3 


2 


N 


f pc (KkIi) = (-l) n K + < +n * +1 8 laka < i p \d+g\kf, > f P (IK) 


3 


3 


Y 


f p (HKk) = S ik f P (IK) 
f x (HKk)=S ik f x (IK)+ 

Sa6 < ic\at a abc\k a > J2 p f p (IK)[(xp\ab) - (ap\xb)] 


3 


3 


N 


fx(Kkli) = 6^ J2 a>b < i a \a^ a ai a \k a > J2 p f p (IK)[(ap\bx) - (ax\bp) 




3 


4 


Y 


f x (HKk) = (-l)< +1 E ai , < i<*\a+ a \k a ><ip\a bfJ \kp > f p (IK)(xb\ap) 


3 


4 


N 


U(Kkli) = (-l)< +1 E a6 < i a \a+ a \k a ><i p \a+p\k p > f P (IK)(bx\ap) 


3 


5 


Y 


f cp (HKk) = {-l) n K +n K +l 8 iB k B < i a \a ca \k a > f p (IK) 


3 


7 


Y 


f pc (HKk) = (-l)<+< + ^5 IafeQ < ip\a c p\k p > f p (IK) 


4 


1 


YN 


Ht Kk = (-ir« + < +1 [E p f P (IK)f P (ki) +E a g a (IK)f a (ki)] 
f cp (UKkf = (-lT* + < + < +1 5 laka < irfa+flkf, > f P (IK) 


4 


2 


Y 


f cp (HKk) = (-l) n « +n K +1 5 igk0 < i a \a+ a \k a > f p (IK) 


4 


3 


Y 


f x (HKk) = (-ir%y ab <i a \a ba \k a ><i,3\a+ \k >yj p (IK)(xb\ap) 


4 


3 


N 


MKkli) = J2 ab < i a \a+Jk a >< ^0+ |*„ > Y,J p {IK){bx\ap) 


A 
4 


A 
4 


v 


f p (HKk) = S ik f p (IK) 
f x (HKk) = S ik f x (IK)+ 
Si n k f3 J2ab < ia\ai a aba\k a > J2 p f p (IK)(xp\ab)+ 
5i ck a Eab < i^ip^lk/} > y /p f p (IK)[(xp\ab) - (ap\xb)] 


4 


4 


N 


f x (KkIi) = 6 iaka y, a>b < ip\at a+ p \k > J2 p f P (IK) [(ap\bx) - (ax\bp) 




4 


6 


Y 


f cp (HKk) = {-l) n * +n K + < +1 5 la k a < if}\a c0 \kg > f p (IK) 


4 


7 


Y 


f cp (HKk) = 5^ < i«\a ca \k a > f p (IK) 



Table 4 continued 



Old case (IK) 


New case (IiKk) 


(Ji) < [Kk)l 


New MEs 


5 


1 


YN 




5 


3 


Y 


f x (HKk) = 5 igkg £ n < i a \a+ a \k a > g cx {IK) 


5 


5 


Y 


f pq (HKk)=S ik f Pq (IK) 


6 


1 


YN 


HttKk=Ha^9ab{IK)f ab {ki) 


6 


4 


Y 


f x {HKk) = (-l) n K +n K +n *S laka < ip\a+ g \k,3 > g cx {IK) 


6 


6 


Y 


f pq (HKk)=6 ik f pq (IK) 


7 


1 


YN 


H£ a Kk =^ nh g ab (IK)f ab (ki) 


7 


3 


Y 


f x (HKk) = {-lY* + < +n * +1 8 laka £ c < i g \&t 8 \k p > g xc {IK) 


7 


4 


Y 


f x (HKk) = (-l) n ^+<S igkg J2„ < i a \a+ a \k a > g cx (IK) 


7 


7 


Y 


f vq {IiKk)=8 lk f pq {IK) 



